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Abstract. We briefly review central concepts concerning nonextensive statistical mechanics, 

based on the nonadditive entropy S q = k — ' (q G 1Z; Si = —k22 i Pi hip;). Among others, 
we focus on possible realizations of the g-generalized Central Limit Theorem, including at the 
edge of chaos of the logistic map, and for quasi-stationary states of many-body long-range- 
interacting Hamiltonian systems. 



1. Introduction 

1.1. Entropy 

The ubiquitous concept of energy is associated with the possibilities for the configurations of 
a mechanical system (e.g., the eigenvalues of the Hamiltonian of a quantum system defined in 
a specific Hilbert space). The concept of entropy emerges in an even larger domain, since it 
can be defined for any system, mechanical or not, which admits a set of probabilities for its 
possible configurations. For instance, if we are dealing with a quantum mechanical system, 
the set of probabilities typically is that corresponding to the eigenvectors of the Hilbert space. 
Epistemologically speaking, entropy is one of the most subtle concepts in physics. Entropy 
and energy together constitute the basis on which statistical mechanics — one of the pillars of 
contemporary physics — is constructed. 

The entropy, initially defined by Clausius for thermodynamics, connects the macroscopic 
and microscopic worlds. Its most elementary form is the logarithmic one, first introduced by 
Boltzmann and refined by Gibbs, von Neumann, Shannon, Jaynes and others. For a finite 
discrete set of probabilities {pi} is given by 

w / w \ 

Sbg = -k^Pilnpi \J2P i = 1 ) > W 
i=i \i=i / 



where BG stands for Boltzmann- Gibbs. The conventional constant k is typically taken to be the 
Boltzmann universal constant for thermostatistical systems, or taken to be unity in information 



theory. For the particular case of equal probabilities, i.e., pi = 1/W, Vi, we have 

S BG = klnW, (2) 

carved on stone in Boltzmann grave in Vienna. Expression (1) enables the construction of a 
remarkably useful physical theory, referred to as Boltzmann- Gibbs statistical mechanics. 

Many entropic forms have been introduced, and reintroduced, since Boltzmann (see [H [2] 
and references therein for details). In 1988 [3], a more general form, namely 

1 _ v w r) q ( w \ 
S i = k g i (5> = 1; « € ft; S x = S BG J , (3) 

was proposed as the basis for generalizing BG statistical mechanics (into a theory now known as 
nonextensive statistical mechanics^). For the particular case of equal probabilities, expression 
d3j becomes 

W 1 ~ q - 1 

S q = k — . (4) 
By introducing the q — logarithmic function 

x 1 ~ q — 1 

ln„ x = (x > 0; q £ 1Z; lni x = In x) , (5) 

1-q 

expressions d3J and Q can be respectively rewritten as follows: 

W w 
S Q = kJ2piln q — = -k^pjlugPi, (6) 
i=i Pi i=l 

and 

S q = k\n q W. (7) 

The axiomatics associated with this entropy have been quite explored and interesting 
characterizations have emerged. For details we may address the reader to [H El El [TJ [8] , among 
others. 

1.2. Additivity versus extensivity 

We adopt for entropy additivity the definition given in Penrose's classical book [9], namely that 
an entropy S is said additive if, for any two probabilistically independent systems A and B, i.e., 
for pfj~ B = pfpf ') V(i, j), we have that 

S{A + B) = S(A) + S(B) , (8) 

where S(A + B) = S{{pf+ B }), S{A) = S({pf}), and S(B) = S({pf }). 

^From definition ([3|) it is straightforward to prove that, for any two probabilistically 
independent systems A and B, 

S Q (A + B) = SJA + ^{B) + S q (A)S q (B) _ 

A. /> n A n 



The word nonextensive is to be associated with the fact that the total energy of long-range-interacting 
mechanical systems is nonextensive, in contrast with the case of short-range-interacting systems, whose total 
energy is extensive in the thermodynamical sense. 



Therefore, Sbg is additive, whereas S q (q ^ 1) is nonadditive 0. 

Entropic extensivity is a concept in some sense more subtle than additivity. An entropy S of 
a given system constituted by JV elements is said to be extensive if 

< lim < oo , (10) 

i.e., if S(N) oc JV for JV >> 1. We see therefore that additivity only depends on the specific 
mathematical connection between the macroscopic entropy functional and the probabilities of 
the configurations of the system. Extensivity depends on this but also on the specific system, 
more precisely on the nature of the correlations of its elements, and therefore of its collective 
configurations. 

The distinction between additivity and extensivity has already been illustrated in simple 
probabilistic systems [10J. It has also been shown for the so-called block entropy (entropy of a 
subsystem of the entire system) of strongly quantum entangled fermionic and bosonic systems 

For the probabilistic system it has been shown [10J that 

S qmt (N)<xN (N»l), (11) 

where ent stands for entropy, JV is the number of (strongly correlated) binary random variables, 
and ^ 

q ent = 1 - - , (12) 
a 

d = 1,2,3,... characterizing the width of an infinitely long strip of nonvanishing probabilities 
of a probability triangle asymptotically satisfying the Leibnitz rule (i.e., asymptotically scale- 
invariant). 

For the (one-dimensional) fermionic system it has been shown that, at criticality at vanishing 
temperature, we have [11] 

lim S gent (N,L)*L (L»l), (13) 

iV— >oo 

where L is the number of first-neighboring spins (or analogous elements) within an infinitely 
long (JV — > oo) chain, and 

qent = , (14) 

C 

where c is the central charge (c = 1/2 for the Ising ferromagnet, and c = 1 for the isotropic XY 
ferromagnet, in the presence of a critical transverse magnetic field in both cases). We verify 
that q en t monotonically increases from zero to one when c increases from zero to infinity. 

For the (two-dimensional) bosonic system the results are qualitatively the same. However 
they have been established only numerically, not analytically. 

The generic scenario which emerges is that, for a vast class of systems (but certainly not all), 
a value q en t exists such that S qent (N) oc JV (JV — > oo), where JV is the number of elements of 
the system under consideration (which might be the entire system, or only a large part of it). 
For standard systems, we have that q en t = 1; for various classes of anomalous systems, we have 
q en t 7^ 1. The situation is depicted in Fig. [TJ 



2 We remark that additivity is obtained whenever (1 — q)/k — » 0. We see that this can occur in two different 
manners: q — ► 1 for fixed k, and k — > oo for fixed q. The latter corresponds to the infinite temperature T limit, 
since in all thermostatistical systems T always appears in the form kT, thus having the dimension of an energy. 
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Figure 1. Entropic additivity and entropic extensivity are different concepts. The BG entropy 
is additive, whereas S q (for q ^ 1) is nonadditive. Extensivity depends on the system: for 
standard systems, Sbg is extensive, whereas it can be nonextensive for anomalous systems; S q 
is the other way around. 



1.3. Central limit theorems 

The Central Limit Theorem within theory of probabilities basically states that the sum of a 
large number N of independent (or quasi-independent in some specific sense) random variables 
whose variance is finite converge, after appropriate centering and rescaling, to a Gaussian (i.e., 
p(x) oc e~@ x ). This distribution constitutes an attractor in the space of distributions, and 
is therefore thought to be the reason for the ubiquity of Gaussian distributions in nature. If 
the single distribution has a divergent variance instead (and also satisfies some supplementary 
mathematical conditions), the attractors are the celebrated Levy distributions. The situation 
changes drastically if strong correlations exist among the N random variables. Depending on 
the nature of the correlations very many types of attractors might emerge. There is however a 
special class of strong correlations, referred to as q-independence [12], for which the attractors 
are g-Gaussians (i.e., p(x) oc e~^ x , where the g-exponential function is defined as the inverse 
of the g-logarithmic one defined in Eq. ([5])), if a specific generalized variance is finite. If this 
variance diverges instead, the attractors are the so-called (q, a)-stable distributions: see [13] for 
full details. The schematic description of these four theorems is presented in Fig. [2j 

The physical-mathematical interpretation of the class of strong correlations named as q- 
independence is not yet fully ellucidated. However, it might well be that g-independence between 
iV random variables implies (strict or asymptotic) probabilistic scale-invariance in the sense that 

J dx N h N (x 1 ,x 2 , ...,x N ) ~ h N -i(x 1 ,x 2 ,...,x N -i) (N -> oo) . (15) 

Although probably necessary, this property is surely not sufficient. Indeed, (strictly or 
asymptotically) scale-invariant probabilistic models (with finite values for the appropriately 



generalized variance) have been analytically solved, some of them yielding, in the N —* oo limit, 
g-Gaussians, whereas other models yield distributions which numerically are amazingly close to 
g-Gaussians, but which definitively are not exactly q-Gaussians. Models that yield g-Gaussians 
are available in [13 [16]; models that have been proved [17] to be not exactly g-Gaussians are 
presented in [181 US] • 

Like Gaussians, g-Gaussians also are ubiquitous H in natural, artificial and social systems. 
What could be the cause of such fact? It could very well be precisely the theorem appearing 
in Fig. [5] which corresponds to g-Gaussian attractors. In what follows we shall exhibit various 
nearly g-Gaussian distributions: in Section 2 for dissipative one-dimensional dissipative maps, in 
Section 3 for long-range-interacting many-body classical Hamiltonians. These two systems share 
a crucial property, namely that they have a maximal Lyapunov exponent which approaches zero, 
thus excluding strong chaos. Finite-size or finite-precision effects are present in them: we mimic 
this property in Section 4 with a simple mathematical model. Finally we conclude in Section 
5 by mentionning various systems presented in the literature which also appear to exhibit q- 
Gaussians. 

2. Unimodal one-dimensional dissipative maps 

Let us here concentrate on a paradigmatic one-dimensional dissipative dynamical system, namely 
the logistic map, defined as xt+i = 1 — axf, where a is the map parameter (0 < a < 2), 
— 1 < %t < 1) and t = 0, 1, 2, .... Our object of interest is the sum 

y 

in the vicinity of chaos threshold a c 
(typically Nq >> 1) and 

(x) 

is the average over a large number of N iterates as well as a large number m n i of randomly 
chosen initial values x± of iterates of the map. This problem has been addressed firstly in [2D] 
and a closer look has been given very recently in [21] . Here, we try to further clarify the study 
and analyse it in a more compact manner. 

Generically, the problem at hand is the form of the probability distribution of the random 
variable given in Eq. (|16p . Indeed, the ordinary Central Limit Theorem (CLT) (yielding the 
Gaussian form) is applicable when the Lyapunov exponent is positive, but, when approaching 
the edge of chaos from above (i.e., for a > a c ), an infinite number of values of a accumulate 
which violate this condition. The situation becomes then quite subtle, as we shall review here. 
Essentially, strong correlations between the iterates of the map emerge. In |21] it is shown that 
the problem is much more complex than the ordinary CLT, which is the case when the map is at 
the chaotic regime (e.g., for a = 2). One needs now to be careful on how we approach the chaos 
threshold point (a c ) and how the number of iterations (N) to be used increases to infinity. In 
mathematical language, this means that two limits are to be performed simultaneously, namely, 
(a — a c ) — ► and 1/N — ► 0. It is argued in |21] that the limit distributions appear to be of 
(/-Gaussian type if these two limits are performed simultaneously in the following special way 
(which satisfies an appropriate scaling relation). We first choose a value of a above and close to 

3 The word ubiquitous is here used not in the strict sense of being everywhere, but only in the loose sense of 
being found very frequently. 
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Figure 2. TV^^^^'l-scaled attractors !F{x) when summing N — > oo ^-independent identical random 
variables with symmetric distribution f{x) with Q- variance ctq = dx x 2 [f(x)]®/ dx[f(x)]® 
(Q = 2q — 1; qx = (1 + g)/(3 — q); q > 1). Top £e/t: The attractor is the Gaussian sharing with /(x) 
the same variance o~\ (standard CLT). Bottom left: The attractor is the a-stable Levy distribution which 
shares with f(x) the same asymptotic behavior, i.e., the coefficient C a (Levy-Gnedenko CLT, or a- 
generalization of the standard CLT). Top right: The attractor is the g-Gaussian which shares with f(x) 
the same (2q — l)-variance, i.e., the coefficient C q (g-generalization of the standard CLT, or g-CLT). 
Bottom right: The attractor is the (q, a)-stable distribution which shares with f(x) the same asymptotic 
behavior, i.e., the coefficient C qa (^-generalization of the Levy-Gnedenko CLT and a-generalization of 
the g-CLT). The case a < 2, for both q = 1 and q ^ 1 (more precisely q > 1), further demands specific 
asymptotics for the attractors to be those indicated; essentially the divergent g-variance must be due to 
fat tails of the power-law class, excepting for possible logarithmic corrections (for the q = 1 case see, for 
instance, [2] and references therein). 



a c . We then calculate the quantity n defined as follows: 

n = : ; , (18 

log 5 

where 5 = 4.6692011... is the Feigenbaum constant. We then denote by k the nearest integer 
value of 2n, and define N* through 

N* = 2 k . (19) 

The g-Gaussians numerically appear to gradually emerge when we choose N = N* and keep 
making (a — a c ) — ► (hence N = N* — > oo). Let us refer to this region as the g-Gaussian 



probability distribution functions (PDFs) one and analyse its borders. To do this, typical values 
for a have been used: they are given in Table 1, as well as their related parameters; a values 
can easily be taken so that 2n values would be obtained with the same precision. Each group 
with the same precision enables the construction of a linear curve in the space 1/N vs (a — a c ) s , 
where s = In 4/ In 5 ~ 0.9: see Fig. 3. Numerical inspection has shown that no other g-Gaussian 
linear curves occur at the left of the largest slope and at the right of the smallest slope in Fig. 
3. All g-Gaussian lines appear to exist only between these two extremes; see examples in Fig [J] 
as well as Fig. [S] where (q,@) pairs of all studied cases are plotted. 

Table 1 . The values of map parameter a and N* used in this work. The values of 2n (obtained 
from the scaling relation), and the values of q and (5 (estimated from simulations) are also listed. 
In this work, the critical value a c is approximated as a c = 1.4011551890920505. 
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The g-Gaussian-like PDFs are lost as soon as the scaling relation (Eq. (fTSj) ) is ignored, i.e., if 
values for N larger or smaller that iV* are used. Two new regions emerge. If we use values for 
N that are sensibly smaller than N* , or, in other words, the value of a that is being used is too 
close to the critical point, peaked PDF's are observed (see also |22j). On the other extreme, if 
we use values for N that are sensibly larger than iV*, we observe in most of the cases Gaussian 




Figure 3. View of the parameter region 1/N vs (a — a c ) s with s = In 4/ In 5 ~ 0.9 coming 
directly from the scaling relation (|18p . Typical values of a are shown, and their correspondiong 
values of N* are determined from the scaling relation (|18p . Any other possible choices for TV* 
yield lines that remain between the lines with the largest and the smallest slopes shown in the 
figure. If one approaches the critical point (origin) along any of these lines in this region, the 
probability distribution function (PDF) appears to gradually approach, excepting for a small 
oscillating contribution, a g-Gaussian. The regions at the left of the largest slope and at the 
right of the smallest slope are not accessible as far as iV* values are concerned. Four of the seven 
g-Gaussian examples are presented in Fig. the peaked (almost vertical line) and the Gaussian 
(almost horizontal line) examples are presented in Fig. 



PDF's, i.e., P(y) = e y " ! /( 2 °" 2 ) /^2ii(j 2 ). Two representative examples for the two regions outside 
the g-Gaussian one are indicated with dashed lines (magenta) in Fig. 3 and illustrated in Fig. 6. 

3. Long-range-interacting many-body classical Hamiltonian systems 

Let us consider the following classical Hamiltonian of N interacting planar classical rotators on 
a (i-dimensional (simple hypercubic) lattice [23] : 

8=1 Z Z i,j 'ij 

where rij runs over all possible distances within the d-dimensional lattice. The particular case 
a = is referred to in the literature as the HMF model [24J, and has been intensively studied 
in the last decade (in its standard representation the coupling constant is divided by N, which 
artificially makes the total energy extensive in the thermodynamical sense). Its dynamical 
molecular approach has exhibited a variety of interesting phenomena: see, for instance, [1] and 




5 10 lu -10 -5 

yP(0) yP(0) 



Figure 4. Data collapse of PDF's for the cases N = N* = 2 k for four of the seven representative 
examples given in Table 1. The (possible) g-Gaussians are approached through finite- N effects 
analogous to those exhibited in Fig. [9j 



references therein. In particular, for the isolated system (microcanonical ensemble) at energy per 
particle equal to 0.69, long-standing quasi-stationary states (QSS) emerge when certain classes 
of initial conditions (usually called water-bag initial conditions) are used. Within the water-bag 
conditions, one may consider initial magnetization equal to zero (usually referred to as M = 0), 
or equal to its maximal value (usually referred to as M = 1), or values in between. Such choices 
influence the specific trajectory of the full system within its 2A r -dimensional phase space (Gibbs 
r space). It has been shown that, for M = 1 initial conditions (possibly for virtually all values 
of initial M), ergodicity is broken. Indeed the summed (over n equidistant instants) one- velocity 
marginal PDF differs when we take ensemble-average or time-average: see [25] and references 
therein. Many of the time-averaged PDF's numerically approach a g-Gaussian (see Fig. [7J for 
one such example). Although no analytical proof is available at the present time, this might be 
a consequence of the g-Central Limit Theorem [12] , within which g-Gaussians are the attractors 
in the space of PDF's. The finite-size effects are illustrated in Fig. [7J 

4. Simple mathematical model for crossover between q-Gaussians 

As a mathematical simple model for finite-size or finite-precision effects exhibiting the gradual 
approach to q-Gaussians, let us consider the following differential equation [T] [26]: 



-JL = - ay r -by* (b > a > 0; q > r; y(0) = 1) 



(21) 



p 




2-q Kn =2- 0.2445... = 1.7555... 



Figure 5. The parameters (q,/3) corresponding to the seven g-Gaussians indicated in Fig. 3. 
These specific seven examples appear to exclude the value 2 — q sen = 1.7555..., which could have 
been a plausible result. At the present numerical precision, even if quite high, it is not possible 
to infer whether the analytical result corresponding to the present observations would be only 
one or a set of g-Gaussians, assuming that exact g-Gaussians are involved, on top of which a 
small oscillating component possibly exists. 



If a = the solution is 

y = e~ bx . (22) 

If b = the solution is 

V = e~ ax . (23) 

If b > a > a crossover occurs from the g-exponential solution for x not too large to the 
r-exponential solution for x large enough; increasing size or increasing precision for specific 
models (such as the logistic map at its edge of chaos, or the HMF Hamiltonian model at its 
QSS state) acts analogously to decreasing values of a towards the limit a = 0, for fixed b. Let 
us address this interesting case. ^From (|2ip we obtain 
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where 2 Fi is the hypergeometric function. In general, this function does not admit an explicit 
expression in the form y{x). An exception is the r = 1 case, which yields 

1 
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(25) 




Figure 6. Representative examples from the peaked region (left) and the Gaussian region 
(right). 
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Figure 7. Velocity distribution for the HMF model at its quasi-stationary state (QSS) [Data 
from Fig. 2 of [25]]. The finite- N effects are visible, gradually approaching a (possible) q- 
Gaussian when both N and n are diverging. An analogous behavior is shown in Fig. [9j For 
comparison, a Gaussian is shown as well. 



The r = case must be handled through the explicit x(y) form, namely 
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(26) 



Let us address now the case of the q-Gaussians. Following the form of Eq. (f2Tj) . we consider 



dy 
d(x 2 



-a r y 



(a q - a r )y q (a q >a r >0;q> r; y(0) = 1) . 



(27) 



If a r = 0, or equivalently if r = q, the solution is given by the g-Gaussian y = e q 



If 



a q , the solution is given by the r-Gaussian y = e T 



For the case a q > a r > and 



q > r, we obtain a crossover between these two solutions, the |x| — > oo asymptotic one being 
the r-Gaussian behavior. 

For r = 1 and q > 1, the solution is given by 
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(28) 



The general behavior of these solutions is given in Fig. [HJ It is evident from this figure that 
this solution is not the most appropriate one for the behavior observed in the neighborhood of 
the logistic map edge of chaos. We notice concomitantly that the appropriate solution for the 
logistic map seems to be very close to the one with r = and q > 1, whose solution is given by 
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with a = ao and b = a q — a®. Indeed, this solution seems to be a very good approximation for 
the behavior of PDFs obtained numerically for the logistic map. This can be seen immediately 
whenever the representative example of the solution given in Fig.[9]is compared to the case given 
in Fig. H^,. 
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Figure 8. 

the case r 
q = 2.35. 



Crossover from ^-Gaussian to successively distant Gaussians, i.e., illustrations of 
= 1 and q > 1. Left: q = 1.7; for comparison, a Gaussian is shown as well. Right: 



5. Final remarks 

Many g-Gaussian-like distributions have been observed in recent years in natural, artificial and 
social systems. Obviously, only for a mathematically formulated model, the hope exists to 
possibly prove analytically that the relevant distribution precisely, and not only approximatively, 




is a g-Gaussian. In all other cases, we can only expect for increasingly high-precision indications 
from real experiments or observations. Computational evidence can and does provide important 
hints, however never a proof. 

This said, let mention in what follows some of the many other systems where g-Gaussians 
have been used to approach the observed PDF's: (i) The velocity distribution of (cells of) Hydra 
viridissima follows a q = 3/2 PDF [2?]; (ii) The velocity distribution of (cells of) Dictyostelium 
discoideum follows a q = 5/3 PDF in the vegetative state and a q = 2 PDF in the starved state 
|28j ; (iii) The velocity distribution in defect turbulence [29J ; (iv) The velocity distribution of cold 
atoms in a dissipative optical lattice |32j ; (v) Velocity distribution during silo drainage [30} I3T] ; 
(vi) The velocity distribution in a driven-dissipative 2D dusty plasma, with q = 1.08 ± 0.01 
and q = 1.05 ± 0.01 at temperatures of 30000 K and 61000 K respectively [33J; (vii) The spatial 
(Monte Carlo) distributions of a trapped 13e Ba + ion cooled by various classical buffer gases at 
300 K [31] ; (viii) The distributions of price returns at the stock exchange [33 [Ml [37] ; (ix) The 
distributions of returns of magnetic field fluctuations in the solar wind plasma as observed in data 
from Voyager 1 [38J and from Voyager 2 [39J ; (x) The distributions of returns of the avalanche 
sizes in the Ehrenfest's dog- flea model [40] : (xi) The distributions of returns of the avalanche 
sizes in the self-organized critical Olami-Feder-Christensen model, as well as in real earthquakes 
|41j ; (xii) The distributions of angles in the HMF model [42J; (xiii) The distribution of stellar 
rotational velocities in the Pleiades [43j . Some indirect evidence is available as well: although 
no (/-Gaussian distribution has been directly observed in some relevant physical quantity, a 
^-exponential relaxation has been seen in various paradigmatic spin-glass substances through 
neutron spin echo experiments [44j . 

Clearly, the simplest hypothesis which would explain the ubiquity of q-Gaussians is the 



validity of the g-Central Limit theorem. The involved random variables would, in such case, 
be expected to be g-independent [12]. The present belief is that (probabilistic) scale-invariance 
is necessary but not sufficient for ^-independence. Further studies are needed to clarify the 
applicability of such ideas to the systems mentioned above, as well as possibly others. 
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